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A kinetic model for the Boltzmann equation is proposed and explored as a practical means to 
investigate the properties of a dilute granular gas. It is shown that ali spatially homogeneous initial 
distributions approach a universal "homogeneous cooling solution" after a few coUisions. The ho- 
mogeneous cooling solution (HCS) is studied in some detail and the exact solution is compared with 
known results for the hard sphere Boltzmann equation. It is shown that ali qualitative features of the 
HCS, inciuding the nature of över population at large velocities, are reproduced semi-quantitatively 
by the kinetic model. It is also shown that ali the transport coefhcients are in excellent agreement 
with those from the Boltzmann equation. Also, the model is specialized to one having a velocity 
independent collision frequency and the resulting HCS and transport coefEcients are compared to 
known results for the Maxwell Model. The potential of the model for the study of more complex 
spatially inhomogeneous states is discussed. 

PACS Numbers: 45.70.-n, 05.20.Dd, 51.10.+y 



I. INTRODUCTION 

Many features of granular gases are captured by an idealized system of smooth, hard spheres with inelastic coUisions 
[1,2]. During the past decade considerable attention has been given to this simple system to understand better the 
mechanisms behind observed qualitative differences between real gases and those comprised of grains. Among the 
most productive theoretical tools for analysis used is that of kinetic theory [3] , and more specifically at low density, 
the Boltzmann equation [4,5]. In recent years important conceptual issues, such as the applicability of fluid dynamical 
equations, have been clarified and quantitative methods have been developed for accurate predictions över a wide 
range of experimental conditions. It is fair to say that the Boltzmann kinetic theory is now the primary theoretical 
method for a quantitative description of granular gases. 

While there are accurate and efRcient numerical algorithms for solving the Boltzmann equation [6] , analytic approx- 
imations are more limited and exact solutions non-existcnt. Such analytic results are of considerable interest because 
they complement numerical solutions with a more penetrating explication of the dominant mechanisms involved in a 
specific application. The mathematical complexity of the Boltzmann collision operatör is the limiting factor in making 
progress, so simpler "kinetic models" have been proposed [7]. This approach has been used with great success for real 
gases with elastic coUisions where several exact solutions far from equilibrium have been obtained and shown to be in 
semi-quantitative agreement with the numerical simulations of the Boltzmann equation [8,9]. Recent applications 
of kinetic models to granular gases have yielded similar interesting exact results [10-12]. The collision operatör for 
a kinetic model is constrained to preserve the most important exact properties of the Boltzmann collision operatör, 
such as a special homogeneous solution and the macroscopic balance equations for mass, momentum, and energy. 
Otherwise, the model is chosen for simplicity and tractability. The objective here is to recall one of the first kinetic 
models proposed for granular gases, the BMD model [13], and to generalize it for both a qualitative and a quantitative 
representation of the Boltzmann equation. It will be referred to as the Gaussian model for reasons that will become 
apparent. 

In the following sections the Gaussian model is defined and applied to the simplest cases of homogeneous states 
and weak spatial perturbations of those states. The motivation for this work is to provide a tool for a subsequent 
more detailed study of spatially inhomogeneous states. For example, recent results suggest that the spectrum of the 
linearized collision operatör for a realistle kinetic model could shed important light on the validity conditions for a 
hydrodynamic description [14,15]. Also, a more practical means to describe boundary value problems is desired for a 
more faithful comparison with experiments. Attention is focused here on spatially homogeneous states for an isolated 
system and on transport coefîicients for small spatial perturbations, as a means to compare and contrast models. A 
more detailed application to inhomogeneous states and shear flow will be presented elsewhere. 
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There is extensive current intcrcst in rclated Maxwell models. Several exact results havc been obtaincd recently 
for the homogeneous state of an isolatcd system using these models [16]. In somc mathcniatical rcspccts thc Maxwell 
models are closer to the hard spherc Boltzmann coUision operatör than thc model studied here. However, its predictions 
(e.g., homogeneous cooling state distribution, transport coefRcients) are quite different from those of the Boltzmann 
equation as recounted below. In contrast, the Gaussian model is structurally simpler but with the capacity to give a 
better representation of known results for the Boltzmann equation. 

The basic results known for the hard sphere Boltzmann equation are summarized in the next section. The ideas 
of kinetic modeling and some existing models are reviewed briefly in section 3, and the Gaussian model is defined in 
section 4. It is an extension and synthesis of two carlicr models, the cUipsoidal statistical (ES) model [17] introduccd 
for gases with elastic coUisions to yield the correct Prandtl number, and the BMD model noted above for inelastic 
collisions. The new model is constructed to describe both elastic and inelastic collisions, retain a realistle velocity 
dependent coUision frequcncy, and to yield the correct Prandtl number for accurate transport coefFicients. 

Also in Section 4, the cxact solution to the Gaussian model is obtained for an arbitrary initial homogeneous state. 
It is shown that this class of solutions approaches a universal HCS solution on a time scalc of thc order of several 
collisions. Hence the HCS is the special state for homogeneous granular gases analogous to the Maxwellian for normal 
gases. The properties of the HCS distribution are shown to be similar to those for the Boltzmann equation at both 
small and large vclocities. In particular similar exponcntial dccay occm-s at largc velocities. The special case of a 
velocity independent coUision frequency is studied as well. In that case thc HCS is the same as that for the BMD 
model. The Maxwell models also have a velocity independent coUision frequency. In this special case Gaussian model 
HCS has properties quite similar to those of thc Maxwell models, inciuding algcbraic dccay for large velocities. 

The Chapman-Enskog solution to the Boltzmann equation for small spatial variations around the local HCS applies 
for the models as well. To Navier-Stokes order this solution is characterized by three transport coefîicients. These 
transport coefficients are comparcd for thc various models and the Boltzmann cquation in section 5. The Gaussian 
model "with velocity dependent coUision frequcncy provides an excellent representation of the hard sphere Boltzmann 
results, över a widc range of inelasticity. 

The HCS is the reference state for linear hydrodynamics. Consequently, knowledge of the exact form of this 
distribution function from the kinetic models provides some additional insight there as well. This is iUustrated in 
Section 6 where the velocity dependence of the hydrodynamic modes [15,14] is calculated. The relationship of these 
modes to the fluxes in Green-Kubo expressions for transport coefîicients is also noted. Finally, the results presented 
here are summarized in the last section and some other interesting applications of the kinetic model are suggested. 



II. HARD SPHERE BOLTZMANN EQUATION 

The system considered is composed of N smooth hard spheres of diameter a in a large volume V. If the density 
is sufHciently small, Na^/V « 1, the one-particle distribution function, f{r,v,t), for the number of particles with 
position r and velocity v at time t is determined from the Boltzmann equation [4,5] 

^ vvV(r,v,f) = J(r,v|/(i)). (1) 



The Boltzmann coUision operatör, J, has the form 

J{r,v\f{t)) = -u{r,t,v)f{r,v,t)+ J dv, J da K{d ■ ^')a-^ f{vy ,t,)f{v,w\,t). (2) 
The first term on the right side represents the loss of particles with velocity v at a rate due to the coUision frequency 

u{r,v,t)= j dvı j daif(â-g)/(r,vı,i). (3) 

The second term of (2) represents the gain of particles with velocity v, where {v', v^} are the " restituting" velocities 
that lead to {v, vı} following a smooth, inelastic hard sphere coUision 

v' = V - ^(1 + a-^){d ■ g)â, v'ı = vı + ^(1 + a-^){a ■ g)â. (4) 

Here, â is a unit vector along the line of their centers, and g = v — vı. The parameter a is the coefficient of restitution, 
O < cü < 1, describing the fractional change in the normal component of the relative velocity (ct • g = — aa • g') and 
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hencc thc inciasticity of coUisions (a = 1 corresponds to elastic collisions. The kernel K{a ■ g) is proportional to the 
flux of particles times the differential eross section and is givcn by 

K{â-g) = a'Q{â-g){â-g), (5) 

whcrc is thc Heaviside step function . 

Thc niost important properties of the coUision operatör are those that result from the microscopic balance equations 
for mass, momentum, and energy in a two partide collision. For the collision rules (4) it follows directly that J has 
the following exact properties 



/ 



dv( V ]j{v,v\f{t))=\ o |, (6) 

im(v-u)V \-lnTC 



where m is the mass, n is the density, T is thc tcmperature, and u is thc macroscopic flow velocity 

n(r,f)u(r,i) \ = I d^r \ v 1 f{r,v,t). (7) 
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n{r,t)T{r,t) I J \ ı^(v-u) 




The two zeros on the right side of (6) correspond to conservation of mass and momentum. The last term results from 
non-conservation of energy and imphes the coohng equation for homogeneous states 

T-^dtT = -C, (8) 

where C, is the "cooUng rate" 

j '^^ j dyrıg' f{r,y,t,)fir,v,,t). (9) 

It is easy to verify from (6) that there is no spatially homogeneous steady state for the isolated system, in contrast 
to gascs "with elastic collisions. Instcad, thcrc cxists a spccial sohıtion, thc homogeneous cooling sohıtion (HCS), which 
is assumed to be approached in a few collision times by ali homogeneous initial conditions. It has a scaling property 
such that the dependence on time occurs only through the temperature 



hcsi^,t)^nv^\t)<j)iv*), v*^v/voit), v„it) - v/2T(t)/m. (10) 

In the following vo{t) will be referred to as the thermal velocity in analogy to a gas with elastic collisions. Sub- 
stitution of (10) into the Boltzmann equation leads to a time independent equation for 4'{v*) that must be solved 
self-consistently with the determination of C, from (9) 



where vq is an average collision frequency 



16 , /tF , , ,^ 



Bccausc of thc scaling property of fhcs both C,* and J* arc independent of time. This problem has bccn studicd in 
detail in recent years and only a few results will be quoted here. For v* <2 the solution to (11) can be expanded in 
Sonine polynomials about a Maxwellian with the result [18] 



with 



c(a) f ^ ^, 15 
1 + -V^ I v*^ - + 



(13) 



32(l-a) (l-2a2) , ^ 

- = 81-17a + 30aMl-a) - ^''^ 
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The subscript B has bccn inchıdcd on csia) here to distinguish the value determined by the hard sphere Bohzmann 
equation from that for the models introduced below. The accuracy of the representation (13) at smaU v to within a 
fcw percent for ali a has been confirmed by Monte Carlo simulation [19,22]. The cooling rate calculated from (13) is 
120] 



<- = ;|(l-»^)(l + |=.w)^ 



For asymptotically large velocities (p{v*) has a qualitatively different behavior [18] 



(15) 



^ (^*) -> Ae-'-l^^' , PBİa) = (16) 
The constant Şb arises in the large velocity limit of the collision frequency, 

t,*{v*)^'^^f3BV*. (17) 

Thus the origin of the cxponcntial dccay is the asymptotic velocity dcpcndcncc of the collision frcqucncy. The ovcr 
population in the tail of the distribution, relative to the Gaussian at small velocities, also has been confirmed for 
V* > 2 for ali a by Monte Carlo simulation [21,22]. These and the transport properties of Section 5 are the primary 
results known for the Boltzmann equation with inelastic coUisions. They are the main features captured by the kinetic 
model proposed here. 



III. AN OVERVIEW OF KINETIC MODELS 



It is remarkable that över a century after Boltzmann wrote his kinetic equation for a low density gas, the content of 
that equation remains masked by its complexity. Certainly, a great deal is known about solutions near the equilibrium 
statc; but the mcchanisms controlling nonlincar transport far from c!quilibrİTinı arc stili poorly undcrstood. Significant 
progress has been made in the past twenty years with the development of Direct Simulation Monte Carlo methods 
(DSMC) by Bird [6]. This numerical tool is exceptionally powerful and provides access now to a wide range of 
nonequlibrium statcs for both clastic and inelastic coUisions. For morc dctailcd analytical insight, kinetic models havc 
provided a parallel powerful tool in rarefied gas dynamics. The objective of this section is to give a brief summary 
of the concept of kinetic models and their extension to inelastic coUisions. Although the discussion is limited to the 
Boltzmann equation, it is noted that the same ideas have been applied as well to its dense fluid generalization, the 
Enskog kinetic equation, for both fluids and sohds [27,28]. 



A. Maxwell Model 



The results quoted in the previous section for the HCS are accurate but not exact. To obtain a more penetrating 
investigation of this and other solutions a simplified version of the Boltzmann equation called the Maxwell model has 
been proposed [23] whereby the kernel K(a ■ g) is replaced by a velocity independent kernel K{r, i).Then (3) implies 
that the collision frequency also is independent of the velocity, ^{r, t) = 47rn(r, t)K{r, t). The resulting model for the 
Boltzmann collision operatör becomes [24] 

J(r,v|/(t)) ^ JM(r,v|/(t)) = -ıy{r,t) 

The collision frequency i'{r, t) is a free parameter of the model, its dependence on space and time is due to a possible 
functional dependence on f{r,v,t). To fix ^(r, f) the cooling rate is calculated directly for the Maxwell model with 
the result 

CM = l{l-a^)ıy{r,t). (19) 

The collision frequency is now chosen to assure that the cooling rate for the model is the same as that for the hard 
sphere Boltzmann equation, (m = C> given to good approximation by (15). This requires the choice 



(18) 
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^(r,i) = ^ (l + â^s^")) ^o^""'*)' (20) 

where uq is given by (12). This completely fixes the Maxwell model. 

The approxinıate Maxwell form for the Boltzmann colUsion operatör does not represent any real kinetics due to 
scattering by a potcntial. It is callcd a Maxwcll model bccausc the property of K being independcnt of the velocity 
follows for scattering by Maxwell molecules interacting via an inverse fourth power law potential. However, the model 
described here retains the coUision rules for inelastic hard spheres, (4), and therefore is a hybrid not corresponding to 
any potential. Stili, it provides an interesting and tractable model for which several exact results have been obtained 
recentiy. 

The HCS has been studied for this Maxwcll model as well. For small velocities (l){v*) again has the form of (13) 
except that the coefHcient c{a) is replaced by [24] 

This is significantly diffcrcnt from the small velocity dependence of the hard sphcrc Boltzmann cquation, suggesting 
that the Maxwell model does not reproduce quantitatively the HCS solution for hard spheres. Furthermore, the 
difference is even qualitative at larger velocities. The exact asymptotic behavior from the Maxwell model is 

(l){v*) ^ Av*-'''-"'l (22) 

Thus there is algebraic decay for the Maxwell model in contrast to the exponential decay for hard spheres. The 
exponent k (a) is the solution to a transcendental equation [25,26] 

' + (1 - »^) p - = (T) ^ + 1 ^ (' - (' - î <' - "'0 ^ ■ 

The behavior of cm (cc) and k (a) is iUustrated in the next section. 

The transport coefRcients associated with Navier-Stokes hydrodynamics also have been calculated for the Maxwell 
model [24] . The agreement with those from the Boltzmann equation for hard spheres is only qualitative (see Section 5). 
While the Maxwell model allows interesting and exact solutions, it does not appear to provide a reliable representation 
of the Boltzmann equation for hard spheres and therefore the results obtained from it must be interpreted with some 
care. 



B. Other Kinetic Models 



The Maxwell model, whilc simpicr tlıan the Boltzmann equation is stili quite complex and even for the HCS the 
exact distribution function has been calculated only in one dimension. Historically, for normal gases, a number 
of simpler kinetic models have been applied with great success. More recentiy, these models have been extended to 
gramılar gascs with a similar success in applications. To explain them generically, it is useful to rewrite the Boltzmann 
equation (1) to make the effects of cooling explicit [29] 

I- + V. /(r, v,i) - ICVv • (V/) = J'[r,v|/(i)], (24) 



dt 



with 



J' (r,v|/(t)) = J(r,v|/(t)) - ^(Vv • (V/) , (25) 
where V = v — u is the velocity relative to the average flow. Then the condition (6) becomes 

J'(r,vi/(t))=( O |. (26) 




In addition, there is a null space for J' {r,v\f{t)) 
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J'(r,v|/o(i))=0. 



(27) 



The conditions (26) and (27) are thc same as thosc for the conservation laws and thc cquilibrium state, respectively, 
for elastic coUisions. More generally, (27) defines the HCS in agreement with (11). These two sets of conditions are 
necessary for the macroscopic balance equations (precursors to hydrodynamics) and the "universal" homogeneous 
statc /o- Thc basic idca of kinctic modcls is to replace the actual Boltzmann colhsion operatör by a simpler structure, 
while preserving the properties (26) and (27). 

There are many ways that a kinetic model can be constructed with these constraints. Perhaps the simplest are 
the BGK nıodel(s) [7] 

J' (r, v|/o(t)) ^ -V (r, t) (/ (r, v, t) - fo (r, v, i)) . (28) 

Clearly, (27) is satisfied and the conditions (26) are imposed by requiring that the relevant moments of / and fo 
should be the same 



/ 



dy\ V I (/(r,v,f)-/o(r,v,t))= I O I . (29) 



This imphes that /o is a functional of / so thc apparcnt simplicity of (28) is misleading. For elastic collisions fo is 
taken to be the local Maxwellian for consistcncy witlı the known equihbrium state. In the casc of inclastic cohisions, 
it would seem appropriate to choosc fo as the HCS distribution from the Boltzmann equation. However, since this 
is not known it is more common to choosc again /o as the local Maxwellian. As described in the previous section, 
this is a rcasonable first approximation to thc HCS if the velocities arc not too large. However, it precludes use of 
the kinetic model to study the HCS itself. The collision frequency (r, t) is a free parameter of the model, usually 
chosen to fit one of the transport coefRcients. On dimensional grounds ly (r, t) (x n (r, t) T^/^ (r, t) and therefore also 
a functional of /. 

The Chapman-Enskog solution to the BGK kinetic equation for inelastic collisions has been obtained to derive the 
associated hydrodynamic equations to Navier-Stokes order (see Section 5) [20]. The dependence of ali transport 
coefHcicnts on thc rcstitıition cocfficicnt a is in good scmi-quantitative agreement with that for the Boltzmann 
equation. Howcver, the model suffers from the same well-known problem for elastic collisions of an incorrect Prandtl 
number rjCp/K where Cp = bks/'irn is the specific heat per unit mass, rj is the shear viscosity and k is the thermal 
conductivity. Since thc BGK model has only one parameter v the absolutc value of either the shear viscosity or the 
thermal conductivity is wrong by a factor of approximately 2/3. This can bc corrected by choosing for fo a more 
general Gaussian, with an additional parameter leading to the ES model for elastic collisions [17] . BGK models of 
this type for granular gases have been discussed recently by Astillero and Santos [30]. 

A related but difîerent kinetic model attempts to represent more directly the gain term of the Boltzmann collision 
operatör [13]. Equation (2) is written as 

J' (r, v|/(t)) ^ -v{vy, t) (/(r, v, t) - g{r, V, < | /)) - ^CVv • (V/) . (30) 

The gain functional ^(r, V, t | /) is now chosen for convenience and simplicity to define the model, but restricted by 
the exact conditions (26) and (27). 

In contrast to the BGK model, the condition (27) now provides an equation that determines a non-trivial HCS 
solution. The simplest choice for g{r,V,t \ f) is again a Maxwellian, but with the temperature modified to account 
for the extra term on the right side of (30). For simplicity, applications of this kinetic model to date have also chosen 
a velocity independent collision frequency. In the limit of elastic collisions it reduces to the BGK model. 

The HCS solution can be obtained exactly for this model and, like the Maxwell model, it has algebraic rather than 
exponential decay at large velocities. The transport coefficients for this second model are of comparable accuracy to 
thosc from thc BGK model and suffer from thc same difficulty of an incorrect Prandtl number. In thc ncxt section 
the kinetic model based on (30) is generalized to inciude a velocity dependent collision frequency and a Gaussian form 
for ^(r, Y,t\ f) that can accommodate the correct Prandtl number. 



IV. gaussian kinetic model 



The main objective of the present work is to propose a synthesis of the BMD and ES models and to extend them 
to inciude a velocity dependent collision £requency. This will be referred to as the Gaussian model. Like the Maxwell 
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model, thc Gaussian model admits exact analysis in many intcrcsting cascs, but it is simpicr and capturcs more 
accurately the qualitative fcaturcs of thc Boltzmann equation. In this section thc model is defined and the initial 
value problem is solved exaetly for spatially homogeneous states. It is shown that ali initial states rapidly approach 
a universal HCS. The HCS is then studied and compared with known results for the Boltzmann equation for hard 
spheres. Finally, it is specialized to the case of a velocity independent collision frequency, for comparison with the 
HCS for the Maxwcll model. 

The model is defined by the choice of a Gaussian for ^(r, v,t \ f) in (30) [31] 

+ V • /(r, V, t) = -u{r, V, t) (/(r, v, t) - g{r, V, t\f)), (31) 



.9(r, V,t I /) = ^(r,t)e^^^^^/('"'*)^^ (32) 
The scalar function A{r,t) and symmetric tensor Bij{r,t) are determined in part by the conditions (7) 

/ M, \ 

Kr,V,%(r,V,i |/)= M2 , (33) 

V M3 - fnTC / 

where Mj are moments of the distribution function, weighted by the collision frequency 

u{r,V,t)f{v,v,t). (34) 





Thc choice of a Gaussian is primarily for convcnicncc and simplicity. Howevcr, it can bc undcrstood also as the result 
from Information theory to determine a function when only the moments in (33) are specified (see Appendix A). It 
follows directly that 

M2(r,f) = Afı(r,t)u(r,t). (35) 
For the special case of constant collision frequency these become 

y4^n(det7rB)"^/^ Ma ^ Mıu; -TrB^—. (36) 

3 m 

so the coefEcients of the Gaussian are related to the density, temperature, and flow velocity. This also illustrates that 

thc elements of B are not fuUy determined by thc moment conditions (33). 

It is useful to divide the matrix B^j into a part proportional to the unit matrix plus a traceless part 

Bij = B6ij + Bij, Bij = {Bij - BSij) , B = -Bkk- (37) 

In the following a tilde above a matrix will be used to denote its traceless part. The special case of (36) shows that 
the trace of Bij is proportional to the temperature and therefore a linear functional of /. It is rcasonablc to c;hoose 
the remaining elements of Bij also to have a linear relationship to /. Furthermore, it is required that this traceless 
part should vanish at /o, the solution to (27), which should be isotropic (to agree with the Boltzmann equation) 

Bii= /'^^^l/=/o(/-/o)- (38) 

Thus, the gain g{r,'V,t \ /) is anisotropic only when evaluated for anisotropic states /.This is an implicit definition 
since /o is a function of Bij. However, the form is such that Bij becomes diagonal when evaluated for / = /o, and 
hence so does B~^. This assures that /o, •when it exists, is isotropic. The trace of Bij is a scalar moment of degree 2 

(inciuding the weight factor ı/) when B^j = 0. Consequently, it is suggestive to take B^j as the traceless part of the 
moment of degree 2 of / — /o- The final form for Bij is then 



Bij = B6ij + dvDij (V) (/ - /o) , Dij (V) = m (vtVj - ^^i,' V') 



(39) 
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where y (a) is an undctcrmincd dimensionless quantity independent of the velocity. The conditions (33) and (39) 
completely determinc thc parameters A, Bij . 

ît remains to choosc the coUision frequency ıy{r,V,t) and the cooHng rate Ç{r,t). in principle, these are specific 
functionals of / in thc Bohzmann equation. Here, they are taken to depend on / only through the temperature and 
density. The coohng rate is chosen to be the same as the Boltzmann result (15) 

C(r, *) = ^ (1 - f 1 + ^CB(a)) Mr, t), (40) 



16^,_^ 2 İ7rT{v,t) 



ıyo{r,t) = -n{r,t)a\ (41) 
5 V m 

Similarly, guidancc for the choice of the velocity dependent colUsion frequency z/(r, V, t) is obtained from that for 
the Boltzmann (see Eq(3)) 

ı^B (r, V, t) = na^ ^ dVı /(r, Vı , t) | V - Vı | . (42) 
For small V this goes to a constant, 

UB - na2n{r, t)V, V = ^ j'^yj' , (43) 
while for large V it becomes linear in V. 

Vb 7r£72n(r, t)V. (44) 

A represcntation of the complctc velocity dependence for the Gaussian model, preserving these limiting forms is 
obtained from (42) using a Maxwellian for / 



v{Y,'V,t) = x{a)vo{Y,t)vlı{v*), vlı{v*) = 



{2v* +v*-^)^eıî{v*) 



(45) 



Here x{a) is a sccond undctcrmincd dimensionless constant. Thc particular choice for x (a) and thc rcsulting accuracy 
of the transport coefBcients is discussed in the next section. It is found that x{a) is a smooth function of a of order 
unity. This form for u*{v*) has the correct large velocity dependence of (16) but with the coefRcient differing by a 
factor of x{a) from thc Boltzmann cquation. 

At this point thc Gaussian model has been spccificd in tcrms of the two remaining constants x{a) and y{a). In 
the next section it will be shown that the thrcc transport coefficients at Navier-Stokes order are functions of two 
independent collision integrals. The constants x{a) and y{a) are chosen to assure that these two coUision integrals 
are the same as those from the Boltzmann equation. This leads to a coupled pair of equations ((110) and (111) below) 
that are solved numerically. This completes the definition of the Gaussian model. 



A. Spatially homogeneous states 

In the rest of this section, attention is restricted to spatially homogeneous states. It is shown that for any arbitrary 
homogeneous initial condition, thc solution goes ovcr to a univcrsal HCS in a fcw collision timcs. For thcsc initial 
conditions u(r, i) = u is constant and by a Galilean transformation it is possible to choose u = 0. Also from the 
continuity equation n{t) = n is constant. The temperature obeys the cooling equation (8) which is now written 

T-'dsT = -C, ds = Uü{t)dt, C = ^- (46) 

The new time variable s represents the average number of coUisions in the time t. It also follows from the definition 
of C(i) that C* is constant. The temperature therefore has a simple exponential dependence on the collision number 

T(s) = e-f**T(0). (47) 
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Now considcr a general hortıogeneous initial distribution and look for solutions to the model kinetic equation in the 
dimensionless form 

fiv,t)=nv^'r{v*,s), v*=v/vo{t), ^* = ^- (48) 
The dimensionless form for the model kinetic equation for homogeneous states becomes 

ds + U* (3 + V* • VvO + ^* {vl] r (v*, s) = u* (v*) g* {v\s\ f*) , (49) 



2 

with 

g* (v*,s I n = «o5(v,i I /)/n = A*(s)e-"' . 
The moment M2 vanishes since u = 0. The remaining dimensionless moments are 

Mı (t) 



Jİ,h/<'v-(.„'..).>w,.)^ (50) 



3nTvo{t) 

The parameters A* (s) and B*- (s) are related to these by 



B*j = B*Sij +yjdv* (v*v* - ls,,v*'^ {f* - /o*) • (52) 

The formal solution to the kinetic equation is found in Appendix B with the result 
r (v*, s) = e-ii'^K{v*,s)rie"^^'^w*,0) 

+ [ ds'e-i^'^' K{v*,s')u*{e-i^'''v*)g*{e-i^'''v*,s- s'), (53) 
Jo 

where K(v* , s) represents the dynamics due to the loss tcrm alone 



K{v*,s)=exp(^- ds"u* (e-5«*«"î;*)^ . 



(54) 



The collision frequency is a monotonically increasing function of the velocity so ly* {v*) > f* (0). This gives the 
inequality 

i^(w*,s) < e-"*W^ (55) 

Since v* (0) is of order unity, the domain of integration in (53) is exponentially bounded for s > 1 and for large s the 
integral becomes independent of s. The first term vanishes exponentially fast and the s independent HCS solution is 
obtained 

f*{v*,s)^(l){v*)= / ds' e- i^' ''K{v*,s')ıy* {e- ^^'''v*)g* {e- ^^^''v*, 00) (56) 
Jo 

It is readily verified that (j^iv*) is a stationary solution to (49) and hence 0(1^*) = /o (f*) is the unique HCS solution. 
It follows from the definition of B that this distribution becomes isotropic on the same time scale as (54) even if the 
initial distribution was not 



B*^{s)^ -Blkicx,)Sij=B*Sij. 



For the class of homogeneous initial states considered, this result shows that the HCS is the universal solution after 
several collisions. Hence it is the special characteristic solution for homogenous states analogous to the Maxwellian 
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for clastic coUisions. The result (56) is strongcr than thc H-theorem for elastic collisions in thc sense that it implics 
the approach to the HCS is pointwise in velocity space. It is interesting to observe that this analysis does not require 
the explicit form for u*{v*)g*{v* , s \ </>*) and so it apphes to models with choices other than thc Gaussian. In fact, it 
applies to the Boltzmann equation itself although in that case (56) is a more imphcit functional relationship whose 
solution must be proved. The consistency of the moment conditions (50) and (51) is verified in Appendix C using the 
exphcit form (56). The functions A* = A* (oo) and B* = B^^. (oo) /3 are fixed by the fact that B*- = O and 

The expUcit forms for these equations also are given in Appendix C. 

B. The Homogeneous Cooling State 

More exphcit properties of the HCS are easily obtained. First, it can be written in the more convenient form 




The coefîicients A* and B* are dctcrmined from (57). Both arc smooth functions of a with thc hmiting values 
A* = 7r~^/^ and B* = 1 ai a = 1. Practical fits for othcr vahıcs of a in thc rangc 0.4 to 1 are given by 

A* = 0.547 - 0.274a - 0.094^2, B*-^ = 1.84 - 0.275a - OMSa"^. 

Similarly, fits for x{a) and y{a) are found from (110) and (111) to be 

x{a) = 0.533 + 0.156a - 0.302a^ 

y{a) = -0.906 + 2.666a - 0.724a^ 
For small v* the form of (l){v*) is Gaussian 

^ö(v*) ^ Cıexp-C2î;*2, (59) 

with 

^ _ A*p _ İ3+P) ( , 3 d^v* \ 

^1-^' - 2{-i+p)v* (0)^^;^'"*=^' 

.= ^^. (61) 

The coefhcient cg (a) in the polynomial expansion (13) for small v* is shown for comparison with the corresponding 
Boltzmann hard sphere result in Figüre 1. A practical fit for cg in the range of a mentioned above is given by 



CG (a) = 0.247 + 0.865a - 2.907a2 + 1.793a^. 



Thc model is seen to rcproduce quite well the Boltzmann results for a > 0.8 and has the same qualitative behavior for 
smallcr a. As indicated in Figüre 1, the corresponding results for both the Maxwcll model and the Gaussian model 
with velocity independent collision frequency (to be discussed in the following section) are always positivc and much 
larger. This is the first of several observations showing an improvement of the model due to the velocity dcpendent 
collision frcqucncy. Figurcs 2 shows thc cxact distribution function rcduccd by thc Maxwcllian at a = 0.9 . Also 
shown are the results from the polynomial expansion using cg (a) and using CB(a) for the Boltzmann equation. It 
is seen that the polynomial expansion follows the exact HCS closely and is close to the expansion using cb{oi). The 
polynomial cxpansion using csia) is vcry closc to thc actual distribution obtained by Monte Carlo simulation of thc 
Boltzmann equation for v* < 2. Therefore the model gives a good representation of the Boltzmann distribution for 
V* < 2. This is found to be true över the whole range of a. 
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I{v*,c) / 2 r dx' \ 
<^(^ ) = ^;^exp --^ / — (x ) , (62) 



The large v* dependence can be obtained as follows. First, rewrite (56) as 

r 



2A* r , , / 2 dx' _ ^ _ . 



= /(oo, c)-— j ^ dxx^ exp - y —i.* (a;') j ^* (a;) e"^ ^ . (63) 

This form holds for any c<v*. The intcgral satisfics thc bound 

/(oo,c) - /(t.*,c) < /oc^iy* (c)e'-^*"'^', (64) 

where Iq is a constant. Since B* is of order unity, this shows that I{v* , c) I{oo, c) for v* > c » 1 with deviations 
of order exp(— î;*^). On this scale of velocities the distribution function becomes 



r ' dx' , , . 



'^(-*)-^^«-P -|^/ (65) 



This is not yet the exponential decay quoted in (16) for the Boltzmann cqııation. Instead, that exponential decay 
requires a stili larger velocity scale due to the slow approach of u* {x') in (65) to its large x' limiting form 

V* {x') ^ f3Gx' + + order e"^" , (66) 

where fİGİx{a) = Pb is the coefficient of the large velocity limit for the collision frequency given in (17). Thus (65) 
behaves as 

<i>{v*) ^ 7(oo, c) exp(-(A/îg(„* - c) + 3) İn t;* + ^(1 - -)). (67) 

The dominant exponential decay is the same as that for the Boltzmann equation, (16), with only the coefficient Pb 
changcd to /^g- In fact, a similar analysis of thc derivation of that rcsult for thc Boltzmann cquation shows that thc 
same intermediate velocity form (67) applies there as well (with (5g Pb)- The cross över to pure exponential decay 
requires very large velocities. In practical terms for a = 0.8 this form holds to within 0.1% for v* > 6 whereas the 
exponcntial decay has thc same accuracy only for much larger velocities. This is illustrated in Figüre 3 for a = 0.8. 
The derivative of the logarithm of the distribution is plotted so that the initial slope for small velocities is near the 
Maxwellian value {—2v), while the asymptotic large velocity value is the constant coefficient of the exponential decay 
( ° ) shown as a dotted line. The intermediate cross över is seen to be governed by the asymptotic form (65) 

v*d.. In cj,{v*) ^ -1 (^3 + Ij.* (^,*)^ . (68) 

Figüre 3 also shows that this more general form persists to very large velocities before the final exponential decay 
is attained. This cross över form is expected to apply for the Boltzmann equation as well and should be taken into 

account in simulation or cxperimental attcmpts to mcasure the över population at large velocities. 
To explore the limiting form for a ^ 1 it is useful to integrate by parts in (58) to get 

= ^.e-«-"---^^'''<teexp(-l £ ^--(/jj {3x'-2B-V)e-'>--'. (69) 

The second term on thc riglıt side vanishes at a = 1 leaving the expected Maxwellian. For a < 1 the second term 
givcs thc cxponcntial decay at large velocities. In order to dominate thc first term it is ncccssary that v* » 1. 
Since the coefficient in the exponential decay is proportional to (1 — a)~^ the relevant domain for overpopulation is 
V* > {1 — a)~^. Clearly for o; ^ 1 this overpopulation becomes physically insignificant. 

In summary, it has been shown that ali of thc qnalitativc fcaturcs of thc HCS for thc hard sphcrc Boltzmann 
equation are reproduced by the Gaussian model. In the next section it is shown that this quality of the model 
extends to transport properties as well. 
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C. Limiting Case: Velocity independent coUision frequency 



To emphasize the effects of the velocity dependence of the collision freguency, it is instructive to consider the same 
Gaussian model witlı a velocity independent collision frequency, u* (v*) 1/*^ (a). Tlıe HCS for the Gaussian model 
then reduces to that of the BMD model. The single parameter of the model, the constant collision frequency, can 
be chosen to fit the shear viscosity or the thermal conductivity. Due to the choice of x{a) made here the model is 
tailored to fit the thermal conductivity. The functional form of u*^ (a) is given by (101) of the next section. The 
HCS solution (56) simplifies to 



(v*) = A* / 
Jo 



(f>{v*)=A* I dse-'e-^'/Pe-''^"'"""^'\ (70) 



where A* ,B*, M^, Mg and p are given by 

A* = {■kB*)-^''^ , B* = fc^, Ml=M^ = v*, p=2v*b/C*- (71) 
A change of variables, t = e~'^^^Pv*'^p/ {p — 2) allows this to be expressed in terms of an incomplete gamma function 



27r3/2 \^ p 

with 



Hvl = ^ r—^ v^-^^^^h -^v*^ , (72) 



7{x,y)= r dte-H^-'. (73) 
Jo 

Interestingiy, the dimensionless distribution function is entirely characterized by the single constant p = 
2ı/*g (a) /(* (a). Its relationship to a is fixed by the choices of cooling rate and collision frequency 

^ ^ 8(l + İ(l-a) + ^CB (g)) 

5(l-a)(l + ^CB(a)) ' ^ ' 

For small velocities the representation (13) applies with c (a) given by 

c{a) - CGi (a) = J^^- (75) 

Figüre 1 shows a comparison of the coefRcient cgi (ce) with that for the Maxwell model given by (21). They are seen to 
be similar for weak dissipation but the Gaussian model grows more rapidly with increasing dissipation. Of course, this 
difference can be eliminated by a different choice of the parameters for the Gaussian model for a closer agreement to 
the Maxwell model rather than the hard sphere Boltzmann equation. The accuracy of this polynomial representation 
is within a few percent for relatively weak dissipation, comparable to that observed for the hard sphere Boltzmann 
equation. It is clear from this figüre that the HCS for models with velocity independent collision frequencies differs 
from that of the Boltzmann equation and the true Gaussian model at small velocities (see Figüre 2). Also, at small 
velocities, the HCS in (70) can be represented as a Gaussian given by 

Tl p(p+3) *2 

â(v*) -J—Ae~(''+^)(^-^)'" . (76) 
' 3+p ^ ' 

The asymptotic behavior for large velocities is obtained from (72) using the limiting form for the incomplete gamma 
function 

This algebraic decay is similar to that of the Maxwell model, and in contrast to the exponential decay for the hard 
sphere Boltzmann equation. This difference is due to neglect of the velocity dependence of the collision frequency in 
both models. Since the exact solution to the Gaussian model is known the crossover from Gaussian to algebraic forms 
can be determined explicitly. Figüre 4 iUustrates this for a = 0.8. The cross över domain occurs for u* ^ 1, increasing 
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slightly with decreasing a, with no special intermediate behavior. Figüre 5 shows a comparison of the exponent for 
the algebraic decay for the Gaussian model, k{a) = 3 + p{a), with the corresponding result for the Maxwell model 
obtained from the solution to (23). Near a = 1 both exponents diverge as (1 — a)~^ but with different coefficients. 

In summary, the simplc Gaussian model with constant coUision frcqucncy capturcs scmi-quantitatively ali of 
the relevant fcaturcs of the Maxwcll model. It has the additional fcature of demonstrating explicitly the solution for 
ali homogeneous states to show the rapid transition to the homogcneous cooling statc, and the detailed features of 
that state. The algebraic decay at large velocities implies divergence of moments of degree greater than some critical 
value for given a. The evolution of such moments for given initial states can be studied in detail to characterize 
the growing över population of large velocities in the HCS. However, the HCS for both the Maxwell model and the 
Gaussian model with constant collision frequency difîer qualitatively from the Boltzmann result for large and small 
velocities. 



V. NAVIER-STOKES HYDRODYNAMICS 

In this section, states with smooth spatial and temporal variations in the density, temperature and flow velocity 
are considered. These are states for which a macroscopic hydrodynamic description is expected to apply. First, the 
results of the Chapman-Enskog method to solve the kinetic equation is recalled. Next, the transport coefficients 
obtained from this solution are identified exactly and in a first Sonine polynomial approximation. Finally, these latter 
expressions are evaluated for the model and compared with the corresponding results for the Boltzmann equation. 

A. Chapman-Enskog solution 

The hydrodynamic equations for spatially inhomogeneous states are obtained from a special solution to the kinetic 
cquation gcncratcd by the Chapman-Enskog method. The method is quitc general and requircs only the propcrties 
(6) for the collision operatör. Since these are preserved in the Maxwell and Gaussian models the results obtained 
earlier for the Boltzmann equation [20] apply for the models as well. The solution is "normal" in the sense that ali 
spacc and time depcndcnce occurs only through the hydrodynamic fields. To first order in the spatial gradients of 
these fields it is found to be 

/ (r, V, t) = /(°) (r, V, t) + /(D (r, V, t) , (78) 

where V = v — u{r,t) is the velocity relative to the flow field. The first term of (78) is the solution to the kinetic 
equation to zeroth order in the spatial gradients 

^C(°^ (r,t) Vv • (v/(°) (r,V,i)) = J[r, v|/(°)(f)]. (79) 

whcrc the superscript on ^^^^ denotes (9) evaluated with Z^"^. Equation (79) is an equation for the velocity dependence 
of /(°^ (r, V, t) which is the same as that for the HCS distribution of the previous sections. The dependence on r, t 
occurs only through the parameters of the HCS. More specifically, J^"' (r, V,f) is the local HCS obtained from (10) 
by replacing the density, temperature, and flow velocity by their actual values in the spatially inhomogeneous state 

/W(r,V,i)=n(r,i)î;o-3(r,t)0(F*), V* = V/vo{r,t), v<y{v,t) = ^2T{v,t) /m. (80) 

The second term on the right side of (78) is proportional to the gradients 

/(ı) (r, V, t) = ^ • V İn T + ö • V İn n + İ [djUi + diUj - V • (81) 

(There are no contributions from the expansion of ( (r. t) to first order as this vanishes for both the Boltzmann case 
and for the models). The functions ^(V|n, u, T), fî(V|n, u, T), and (V|n, u. T) are solutions to the integral 
equations 

(-C('^)TdT + C-Ç-^A = A, (82) 
(^-Ç^TdT + r) B = B + C^°U, (83) 
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(-C(°)TÖT + r)Ci,=Q,-, (84) 

with the definitions 

A(V|n,u,r)= Q + iv-Vv)/(o)V-V/(o)-^Vy/(°), (85) 

B(V|n,u,T) = -V/(o) - -V^/(o), (86) 

m 

a,(V|n,u,T) = l/i(öy,/W). (87) 
The linear operatör C is the collision operatör expanded to first order in /^^^ 

= /dv- y|^;;^(;j] l/-/<°^ (88) 

B. Transport coefficients 

The Boltzmann cquation and ali models considered give the same macroscopic balance cquations for mass, energy, 
and momcntum (or dcnsity, tcmpcrature, and flow velocity) because thcy ali imply the propcrties (6). The Navier- 
Stokcs hydrodynamic cquations follow by evaluating the momentum flux Pij and the heat flux q in the macroscopic 
balance equations using the Chapman-Enskog solution to first order in the spatial gradients, with the results 



Pij = -r] i djUi + diUj - -öijV • u ) , q = -kSİT - ııVn. (89) 



The first of these is Newton's viscosity law, where rj is the shear viscosity. The second is a generalization of Fourier's 
law, where k is the thermal conductivity. There is an additional contribution for granular gases proportional to the 
dcnsity gradient, with a transport coefBcient /x, that does not occur for normal gases. These are identified from the 
Chapman-Enskog solution as [20] 



(90) 



«=^(î^«-2C(°))-i(l+c), (91) 



with the definitions 

_ /dVAj(V)£C»j(V) _ /cA^S(V)-r.4(V) _ /dVS(V)-£B(V) 

" /rfVA,(v)Cy(v) ' " /rfvs(v)-^(v) ' " /dvs(v)-B(v) ■ ^ ' 

The functions Z)jj(V) and S(V) are 



Dij(V)=m[V^j~-VH,A, S(V)=V(-mV"--r). (94) 



İF'^.,). s(v).v(i,„v--5: 

Also, Ko = 15rio/Am and r/o = 5 (mT)^^^ /ISct^tt^/^ are the low density values of the thermal conductivity and the 

shear viscosity in the elastic limit, rcspcctivciy. The constant c(a) is the same as that occurring in the representation 
(13), appropriate for either the Boltzmann equation or the model being considered. The forms (90) - (92) provide the 
exact expressions for these transport coefücients. 
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C. Sonine Polynomial Approximation 



More explicit results reguire determination of f^^^ and the solutions A, B, and Cij to the linear integral eguations 

(82)-(84). The Gaussian model allows explicit constnıction of these. Iîowever, in general it is useful to represent 
these quantities as an expansion in a complete set of polynomials and generate approxinıations by truncating the 
expansion. In practice the leading term in these expansions provides a very accurate description över the full range of 
dissipation and density. The determination of /'■^•' to leading order in the Sonine polynomial has already been given 
by (13). Similarly, the leading contributions to the cxpansions of A, B, and Cij are found to be [20] 

AV) \ ( c^s(v) 

B{V) \^Îm{V)\ ceS(V) , fM{V) = n{nvl) (95) 



CijiV) \ccDij{\) 



with the coefBcients 



- 15nT3 J [ B{V) ■ S(V) ) " ^ -flM ) 



j . , ( ] = [ ) , (96) 

= ^ / = ~ î''- ^^^^ 

The distribution function /'^^ in this approximation is obtained from ( 81) 



M 



2m 11/ 2 \ 

^ (kVT + MVn) • S(V) + ^77- [diUj + djUi - -5,,V • uj A,(V) 



(98) 



To evaluate the transport coefHcients the forms (90)-(92) are used, with the frequencies v-nioı)-, i'fc(o;), and v^ia) 
determined in this approximation by 

_ JdVDijCfMDij JdVS.jCfMS 

JdVfMDijDij ' JrfV/MS-S • ^^^^ 

These integrals have been calculated for the Boltzmann equation [20] 

= = (ı + 1 (1 - o) + m) ^ (101) 

The average local frequency i'o(r, i) is given by (41). The corresponding results for the Maxwell model [24] are 

(1 + a) (4 - a) 3 



'rıM f, ■ 32 



(1 + ^cs(a)), (102) 



Km = Km = ^ (1 + a) (19 - İla) (1 + ^CB(a)). (103) 

Finally, it is straightforward to perform the same calculations for the Gaussian model considered here. The form 
of the linearized coUision operatör L for the Gaussian model is obtained in Appendix D 

Here P is a projection operatör defined by 

VX = 1/5(0)^. y rfv^aX, ffW = 5(1 /(°)), (105) 
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and {tpa} is the orthonormal set 




3_V 

N2 



^3 V ^1/7 



(106) 



The normalization constants are given in Appendix D. With these results the frequencies i^^c ı^nG^ a.nd t'^c are 
found to be 



(107) 



x/dvs-(ı-7^)i^;;^/MS 

JdVfMS^ 



(108) 



The constants a; (a) and y (a) are now chosen to assure accurate transport coefficients. This is most directly done by 
requiring that the above frequencies are the same as those from the Boltzmann equation, i. e. 



VkG = I^kB- 



(109) 



It follows form (90) and (91) that the Prandtl number at a = 1 is Vy^jvy^. So this choice assumes that the Gaussian 
model also will have the correct Prandtl number in the elastic limit. This gives 



X (") = (a) 



JdVfMS^ 



/dVS-(l-P)î/î,/MS' 



(110) 



2/ (a) = -2 ( ı^nB (a) - x{a) ■ 



JdVfMDijDij ) l JdVg(.o)DijDij 



(111) 



With these choices, the transport coefBcients are given by (90)-(92), and the only differences from the Boltzmann 
values results from the replacement of cb {a) by cg (a) in the expressions for k and /i. It should be noted that Eqs. 
(111) and (110) are implicit since the right sides depend on x{a) through the collision frequency in Eq. (57) that 
dctcrmincs the paramctcrs A and Bij. In practicc the calculation of x{a) is done itcratively. First, the intcgrals 
in (110) are evaluated at a = 1 to determine a zeroth order estimate for a; (a). Then (57) is used to get a first 
approximation to A and Bij. Next, these results are used in (111) and (110) to calculate the first approximation to 
X (a) and y (a). The process is rcpcatcd starting with this first approximation for x {a) . The results reported here 
are for two iterations, showing good convergence of the process. The fits obtained for x and y are 

x{a) = 0.533 + 0.156a - 0.302a2, 



y{a) = 0.906 - 2.666o! + 0.724»^ 



For a = 1 these results reduce to 



Kg (1) = (1) = 1, <G (1) = '^;g (1) = <B (1) = ^, 



(112) 



.(l) = il^, ,(1) = -1^. (113) 

These results define a new kinetic model for normal gases, extending the ES model to one with a more realistle velocity 
dependent collision frequency. 

Also, for the special case of a constant collision frequency the general results reduce to 

x{a)vl, = v:^{a), 2,(a) = -2(^Ç^-l). (114) 
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This casc is relevant also if onc wantcd to usc thc Gaiıssian kinctic model to represent thc Maxwcll model. Tlıen 
in (114) f*^ (a) and u*g (a) should be replaced by (a) and v*j^ (a), respectively. Finally, for both a = 1 and 
constant coUision frequency the usual Ellipsoidal Statistical model is recovered 

'-;g(i) = i, ^:g(i) = '^;g(i) = ^'^m = | y = -ı- (ns) 

Figurcs 6,7 and 8 show thc shcar viscosity, thcrmal conductivity and the ncw transport cocfRcicnt /i for thc various 
modcIs comparcd with thc Boltzmann cquation rcsults. It is sccn that the shear viscosity for thc Gaussian model 
with eithcr a vclocity depcndcnt or velocity indcpcndcnt coUision frequency is indistinguishablc from thc Boltzmann 
result. The small differenccs bctwecn the velocity depcndcnt coUision frcquency Gaussian model and the Boltzmann 
results for the k and jj, coefhcients are due to the differenccs between cb(q) and CQ{a), and by truncation of the above 
iteration solution for x {a) after two steps. The difîerences in the case of the constant coUision frequency models 
are more pronounced because as seen in Figüre 1, the CGi{ct) and CM{ct) are significantly different from CB(a) for 
smaller a values. These results show that the Gaussian model has the ability to fit the transport propcrties to the 
hard sphcrc Boltzmann rcsults for ali a, inciuding thc corrcct Prandtl numbcr rjCp/K = 2/3 at o: = 1. The other 
model do not have this capacity and the associated transport coefRcients do not represent as well those from the 
Boltzmann equation [24], although they yield the correct Prandtl number at a = 1. Clearly, the inciusion of the 
velocity dependent coUision frequency in the model allows excellent agreement with the Boltzmann results. 

VI. HYDRODYNAMIC MODES AND GREEN-KUBO EXPRESSIONS 

The simplest solutions to the Navier-Stokes equations are those for a large system with small perturbations about the 
HCS (not the local HCS as considcrcd abovc). The resulting five indcpcndcnt solutions are rcfcrred to as hydrodynamic 
modes. For a gas with clastic coUisions, these would correspond to shear diffusion, heat diffusion, and dampcd sound 
propagation. The hydrodynamic modes are more complicated for inelastic coUisions but their properties have been 
worked out and discussed [20]. 

The Chapman-Enskog method provides a "normal" solution that implicitly presumes the existence of a hydrody- 
namic dcscription. A more fundamcntal study of the contcxt or validity of hydrodynamics is possiblc by dctcrmining 
the possible solutions to the Boltzmann equation for small perturbations of the HCS. The resulting linearized Boltz- 
mann equation is obtained by substituting / = fh^s [1 + A] into (1) and retaining terms linear in A 

(ö, +v*-V*+/:o)A = 0. (116) 

The dimensionless units of Section 4 have been used, and thc linear operatör Co is defined by 

£oA = r'.C(<^A)-r'|^-(<^A). (117) 

For elastic coUisions the second term of (117) vanishes, </> becomes the Maxwellian, and Co is the usual linearized 
Boltzmann coUision operatör. Its spcctnım inciudcs a fivc fold degcnerate valuc at zcro. The corresponding cigcn- 
functions are A — » linear combinations of l,v*,v*^. These eigenfunctions are known as the summational invariants 
because their sum for two particles is conserved in a two partide coUision. Thus, for a = 1 the eigenfunctions and 
eigenvalues of Ca are known and constitutc thc hydrodynamic modes in thc long wavelcngth limit. 

The Identification of the linear combinations of these hydrodynamic modes as eigenvalues and eigenfunctions of Co 
has been given recently [14,15] with the results 

j0.oXn = KXn, (118) 

Aı=0, A2 = Ş, A3 = A4 = A5 = -^. (119) 

The degeneracy for elastic coUisions is partially broken, with some zero eigenvalues going to ±^o/2. The corresponding 
eigenfunctions are 

Xı=4 + vd^lncP{v*), X2 = -^-vd^\n<l){v*), (120) 
Xn = v„d,ln<j>{v*), n = 3,4,5. (121) 
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For a = 1, dv\n(j> {v*) = —2v and the Xn become linear combinations of l,v*,ı;*^. This suggcsts that (119)-(121) 
provide thc hydrodynamic modcs for a < 1 as well. This is confirmed by noting that thcsc cigcnvalucs arc the same 
as those of thc macroscopic balance equations in the long wavelength Hmit. 

The velocity dependence of the hydrodynamic modes for a < 1 is comphcated due to their definitions in terms 
of the HCS distribution. An advantage of the Gaussian models is that this distribution is known exphcitly and the 
construction of the eigenfunctions is straightforward. AH of thcsc modcs are characterized by dy İn cj). This has aheady 
been shown in Figüre 3 for a = 0.8. Figüre 9 shows the same data but with the result for the velocity independent 
collision frccjucncy model inciuded. The dashed curve in each case represents the elastic a = 1 limit. For the velocity 
independent coUision frcqucncy xı approachcs a constant for largc v, duc to thc asymptotic algcbraic decay of (p. For 
the velocity dependent collision frequency it approaches v according to the exponential decay of (j). 

Figüre 9 shows that there are significant qualitative difîerences from the hydrodynamic modes from the elastic limit 
when V* > 2. This is thc crossovcr of the distribution function to its largc velocity form (68). Thc fact that thc 
hydrodynamic modcs arc rclated to the log of the distribution function lends ncw importance to these asymptotic 
forms. 

Related properties are the fluxes appearing in the Green-Kubo expressions for the transport coefhcients. The 
expressions in the previous section can be written in a form suggestive of Green-Kubo relations for normal fluids [34] 

.= ^^/^rf^'(^^-^U(^'))e-^*% (122) 
n=^^ rds'{S*.^l{s'))e'^i'^', (123) 



f^=l^ + £ ds' (S* • (*î(5') - *^(a'))) . (124) 

The brackets denote an average över the HCS in the dimensionless velocities 

{X) = j dvXv*)X(v*). (125) 

Furthermore, the dependence on the dimensionless time s is defined by 

X{s) = e"'^«X(V*). (126) 

The averages in these expressions therefore have the interpretation of time correlation functions. The momentum 
flux _Dj* and heat flux S* are the same as in (94). They arc fluxcs in the usual sense of the velocity v times linear 
combinations of the summational invariants l,v*,t;*^. In the elastic limit the other functions also have these 
forms 

*î^0, <Î>U-^İ5*, *;^S*(V). (127) 

The resulting expressions (122) and (123) for rj and n are then precisely the low density limits of the usual Green-Kubo 
expressions as time integrals of flux autocorrelation functions [32]. 

For a < 1 thc functions $* are no longer simply rclated to fluxes of the summational invariants. Instead they can 
be written as fluxes for the hydrodynamic modcs dcfincd above 

*î =v*(xı+X2) + ^X, (128) 
*2,y = \ (viXo - \5io^-Xj , i = 3, 4, 5, (129) 

*3 = ^(v*X2 + x), (130) 

where x is the vector whose components are Xn, n = 3,4, 5. This relationship of the "fluxes" to the hydrodynamic 

modcs is thc same as for a normal gas. Only thc forms of thc hydrodynamic modes change for a < l. Howcvcr, since 
these modes are significantly difîerent at large velocities, it is expected that their efîect on the transport coefhcients 
may be important. 
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VII. DISCUSSION 



The objective here has been to describe a simple but reaUstic kinetic model for the hard sphere Boltzmann equa- 

tion. The ncw fcaturcs of thc Gaussian kinetic model defincd in section 4 rclative to previous modcIs are 1) a vclocity 
dependent coUision frequency, 2) two free parameters for a good description of transport coefîicients, and 3) appli- 
cability to both elastic and inelastic collisions. For elastic collisions and constant collision frequency it reduces to 
the ES kinetic model [17,7], while for inelastic collisions and symmetric Gaussian it reduces to the BMD model [13]. 
For elastic collisions, constant collision frequency, and symmetric Gaussian it becomes the usual BGK model [7]. It 
is also shown here that the Gaussian model for constant collision frequency can be "tuned" to represent well the 
more complicated Maxwell models. One motivation for the generalization of a kinetic model to inciude a velocity 
dependent collision frequency is a more accurate description of the overpopulation at large velocities for granular 
gases. The decrease of thc distribution function for large velocities in the simplest statc of HCS is algcbraic for any 
model with a constant collision frequency, inciuding the Maxwell model. In contrast, the decay found from the hard 
sphere Boltzmann equation is exponential due to the velocity dependence of the loss term in the collision operatör. 
This qualitative diffcrence may be important for driven states as well. Although this asymptotic behavior occurs only 
for extremely large velocities it can have an effect on the moments of the distribution function. In addition it has 
been shown in Section 6 that the hydrodynamic modes and the Green-Kubo fluxes depend on the log of the HCS 
distribution function, so this asymptotic behavior is even more important. The Gaussian model with velocity depen- 
dent collision frequency incorporates this behavior and in addition gives a quite good quantitative representation of 
the HCS distribution function for small velocities as well. This is iUustrated in Fig 1 where CG{a) shows significant 
improvement över the velocity independent case. As a consequence the transport coefHcients «(a) and jjbia) are also 
significantly improved due to their dependence on cgİo)- 

The second feature of a non-symınetric Gaussian provides an additional parameter beyond the collision frequency 
that can be chosen to optimize the quality of ali transport coefficients. Here they are chosen such that the shear 
viscosity is accurate for both the constant collision frequency and the velocity dependent collision frequency for ali 
values of the restitution coeflicicnt. The other transport coefficients are accurate in thc elastic limit, inciuding thc 
correct Prandtl number for both cases. For inelastic collisions the agreement with Boltzmann remains excellent for 
the velocity dependent collision frequency case for ali a. This is a primary improvement of the Gaussian model. 
In contrast, the transport coefRcients from the Maxwell model are quite different from those of the hard sphere 
Boltzmann equation, and the other kinetic models using a symmetric Gaussian ali give the wrong Prandtl number. 

An advantage of most kinetic models is their structural simplicity. They can be solved exactly for many states as 
functionals of a few moments of the distribution. These moments stili obey complicated nonlinear integral equations 
but the problem is simplified to the extent that exact results are often possible for states with sufScient symmetry. 
An example is given here for homogeneous states whcre the exact solution is obtained in terms of the parameters of 
the Gaussian gain term, A (s) and (s) , which in turn are defined in terms of the moments M), {s). It is shown 
that an arbitrary homogeneous initial condition evolves after a few collisions to a universal scaling solution, the HCS. 
Such behavior is expected also from the hard sphere Boltzmann equation but its complexity has precluded a proof 
to date. It is useful also to have the explicit representation of the HCS for other purposes as well. Here it has 
been noted that the hydrodynamic modes for weakly inhomogeneous states are described by eigenfunctions of the 
linearized Boltzmann collision operatör. These eigenfımctions are determined from the HCS and a contrast with the 
corresponding eigenfunctions for elastic collisions was made possible by the explicit results for the HCS for the kinetic 
model. Significant differences are observed between the cases of the velocity independent and velocity dependent 
collision frequency, due to the qualitative differences in the large velocity dependences of the HCS. This is also related 
to the Hilbert space for formulating the cigenvalue problem for the linearized kinetic equation. The natural scalar 
product is an integration över the velocities weighted by the HCS distribution function. Due to the algcbraic decay at 
large velocities for the constant collision frequency case (inciuding the Maxwell model), polynomials of high degree do 
not exist in this space [14]. This restriction does not occur for the hard sphere Boltzmann equation or the Gaussian 
model with velocity dependent collision frequency. It is of somc interest to study any qualitative differences in the 
spectrum of the linearized collision operatör and any consequences for the existence of hydrodynamics. The Gaussian 
kinetic model provides a tractable context to address this issue. 

The most interesting states for cxperimental purposes are quasi-steady states for systems driven at the boundaries. 
For states of high spatial symmetry the kinetic model again offcrs the advantage of an exact solution as a functional 
of low degree moments. An example is that of unifornı shear flow where an exact solution for the distribution function 
has been obtained in the case of a symmetric Gaussian [35]. The result applies even for large shear rates so the 
rheology of states far from equilibrium can be studied directiy. The Gaussian model described here also can be 
solved cxactly for unifornı shear flow and will be given elsc;wherc. Vibrated systems, with and without gravity, have 
been studied on the basis of the Boltzmann equation using Monte Carlo simulation methods leading to a number 
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of important rcsults bearing on experimcnts (c.g., boundary laycrs [36], dcpcndencc of vclocity distribution on the 
distance from the driving wall [38], synımctry brcaking [37]). The Gaussian kinetic model may be simplc cnough for 
a complementary analytical study of such problcms. 

In summary, the work here has extended earher kinetic models to bring closer correspondence with the Boltzmann 
equation for the HCS and small spatial perturbations of that state. The price for these improvements is an increased 
complexity of the model, although this has not becn an impcdimcnt for the simple states considered here. It remains 
to demonstrate significant new results for more complex states, not already addressed by the simpler existing models. 
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APPENDIX A: MOTIVATION FOR GAUSSIAN MODEL 



The Gaussian model results from an approximation to the gain contribution to the Boltzmann coUision operatör, 
denoted by .g(r, V, t | /) in (30). The specific choice of a Gaussian can be interpreted as resulting from maximizing 
the Information entropy I[g] 

I[g] =Jdıı (v) g{r, V, i | /) lng{r, V, 1 1 /), d/ı (v) = dvu{v), (Al) 

among the class of functions whose weighted moments of degree 2 are specified 



d/z(v) V \g{r,V,t\f)= G2 . (A2) 

The measure for the velocity integration has been cliosen to incinde the velocity dependence of the collision frequency. 
Incorporating these constraints with Lagrange multiplyers and minimizing I[g] leads directly to the Gaussian form 

g{r, V, i I /) = exp (-Aı - A2 • \-XsijViVj) , (A3) 

where the coefhcients Aa are determined in terms of Ga from (A2). Thus, if the only known or important exact 
properties of the gain term are the moments in (A2) then (A3) is a " natural" choice for the model. 

It may be useful to recall that ^(r, V, i|/) is exactly Gaussian for / = Maxwellian at a = 1. It has been verified 
numerically that this property remains tme to an excellent approximation for a < 1 as well, with only the parameters 
of the Gaussian changing. This gives further support for the choice (A3). 

APPENDIX B: FORMAL SOLUTION FOR HOMOGENEOUS STATES 
The formal solution to the Gaussian model kinetic equation (49) is 

/*(t;*,s) = e-(5^*(3+^*-^-'^+^*(''*')V*(î^*,0) 

+ rds'e-(^«*(^+^*-^-*'+"*(''*))(^-^')ı.*(î;*)5*(î;*,s'). (Bl) 
Jo 

The action of the exponential in (Bl) can be determined as follows. Define a function X{v*,s) by 

X{v\ s) = e-(K*(3+v*-v,.)+'^*(''*))«X (^;*) ^ (B2) 

which then obeys the equation 
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ds + \C (3 + V* • VvO + ı^*{v*)^ X = 0. (B3) 
Next introduce 

X{v*, s) = e-i^'''''-^-'X{v*,s), (B4) 

so that X{v*, s) obeys the equation 

(ds + + e^f*"^'-^"* ı/*(î;*)e-^«*''^*-^-*^ X = 0. (B5) 

Prom the identity 

gir«v*-v^. j^(^*-)g-iC*«v*-Vv. ^ p (^e^'^'-'v*^ , (B6) 

this equation becomes 

(ds + ^C + t^* {ei^''v*yjX = 0. (B7) 
This can be integrated directly and inserted in (B4) to give 

X {v*,s) = e-i^*^e-5f*«^*-^v. exp ds'v* (e^^'^'u*)^ X {v*) 

= e-i^'' exp ds'v* (e-5f*«'î;*)^ X (e-^f* . (B8) 

The fornıal solution to the kinctic cquation bccomcs 

f*{v*,s)=e-i^''K{v*,s)f*{e-^^'''v*,Ö) 

+ J\s'e-i^'''K{v*,s')u* (e-^^'v*) g* (e^^'î;*, s - s') , 



(B9) 



r ds'v* [e-i^'^'v*^^ . 



K{v*,s) ^c^pi- ds'v* le-^^'' V*)) . (BlO) 



It is intcresting to nete that no use of the exphcit form for g* has been used. So, this result appUes to the Boltzmann 
equation as well. 

APPENDIX C: MOMENT CONDITIONS 

The HCS for the Gaussian model is given by (56) 



/•OO 

= / ds'e-i'^'''K{v*,s')v*{e-^'^''''v*)g*{e-^'^''''v*,oo). 



(Cl) 



This is restricted by the moment conditions (50) 



M* 
Mİ 



*)<Piv*). (C2) 



These conditions can be verified by direct integration 



M- 
M3 



İ^= J d^r* ) v*{v*) J^^ ds'e-ii'''K{v*,s')v*{e-ii'''v*)g*{e-i'^'''v*,oo) 
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/ ( |İ*2 )^J^ dxx^u*{x)g*{x,<^)^y{v*)K{v\x), 



with the notation 



Next eliminate 2u*{v*)/C,* by noting it can be generated by differentiating K 



, .dK{v*,x) 
{x,co)- 



= -Att r dxx''u*{x)g*{x,oo) H dv* (2^., ^ ^^(^*'^) 



4w dxx'^ıy*{x)g*{x,oo)i^(^2l.2^ + J dv* i'^^^ K{v* ,x) 

where use has been made of (51). The second term of (C3) can be recognized as moments of (l>{v*) to give the desired 
result 



m; \ _ / Mı* 



Mi J \m;-c 

_ ( M* 
M3* 



(C4) 



The last equahty foUows from (7) in the form 



This confirms the consistency of moment conditions (50) and (51). 

The two equations (C5) fix the values of A* and B*~^ in the Gaussian model (49). A convenient representation is 

3 _ d.v*v* /ff dxx'^v*{x)e-^"^''^K{v*,x) 
2 jo" dv*v*^'^ ^Q dxx'^i'*{x)e~^*~^^^ K{v* ,x) 



01 



0=1 dv* { V* - ^v*-^ ] I dxx^u*{x)e-^' '^'k{v*,x). (C6) 



o \ 2 / jQ 

This determines B*. Next A* is obtaincd from 



A*-^ = ^J dv*v*-' I dxx'^ıy*{x)e-^' '='''k{v*,x). (C7) 



FinaUy, with A* and B* ^ known, the moments are determined from (51) 

Mı* 
M3*-C 



^1 ^ - I dv* ( 2 ^2 ) ıy*{v*)A*e-''' (C8) 
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APPENDIX D: LINEARIZED COLLISION OPERATÖR 



The collision operatör for the Gaussian model is 

J ([r, v|/(0) = -v{v, V, t) [/ (r, V, t) - g{r, V, t \ f)] . 
The distribution function is expanded as 

/ (r, V, t) = /(°) (r, V, t) + /(ı) (r, V, t) + .., 

g{r,V,t\ f) = g{v,Y,t \ /(°)) + j dv'^-^^^^^ |,=,(o) f^'^ (r, V',t) + ... 
The Unearized colhsion operatör is therefore 



Sf{T,v',t) ^ ^ 



u{T,V,t) 



(Dİ) 

(D2) 
(D3) 



(D4) 



The second term of (D3) can be made more exphcit by recalhng that the functional dependence of g{r, V, f | /) occurs 
only through A and Bij 



I 



i) 



(D5) 



Using the definition of Bij in (39) gives 



-B-Hr,t\f^°^)'-^^^B-Hr,t\f^^^) 

9 f^^B,.ir.t\f)^^^^yMn^^^v)). 



(D6) 



Then (D5) becomes 

,'^.9(r,V,t I /) 



/ 



öf{r,v',t) 



(r, V',*) = ,(r, V,t I /(O)) [/ ^v'^^^^^^ |,.,.) /« (r, V',t) 



/ ^v'f^ 1,^.0, /« (r, V',) . J^A.(V) / .v'.,(V')/« (. V',) 



(D7) 



The expansion (D2) leads to an corresponding expansion for the moments 

ir,v,t) (/(O) (r,V,t) + /« (r,V,i) + ..) 



M,\ 1 
M2 = / dv V I jy( 

M3 y V im (v - u)2 



M 

= I M 



(0) 

1 

(0) 



(D8) 



M. 



(1) 



The coefficients A{r,t \ Z^"^) and Bij{r,t \ /") = (Bkk{r,t \ /°)/3) = Bdij are determined from Mj^ just as is 
done in Appendix B. The remaining part of the moment conditions are 
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M 



(1) 



(D9) 



The terms on the right side of (D7) can be identified as an expansion in terms of polynomials of degree 2 in the 
velocity. To do so, first define a Hilbert space with scalar product 



{a,b)= jdW^^a*h, fl(°)=ö(l/^°^) 
Next, define the set of functions {V'o-} 



(DIO) 



with normaHzation constants 

Arı = (l,l), N2 = {Vi,Vi) = {l,V^) 
These functions form an orthonormal set 

Equation (D7) may now be written in the form 



N2 



\ 



1 I y2 N2 



(Dil) 



N2 
İVı 



N2 



Ar2 



dv 



,Sgir,\,t\f) 



/W (r,V',t)=5(r,V,t 1/(0)) 



ea^a (V) + 



nm?B'^ 



The coefHcients can be determined by taking the scalar product of this equation with 



A,(V) j dVD,,f('^ 



(D12) 



(D13) 



(D14) 



(D15) 



The sccond cquahty follows from (D9) and allows these terms to be represented as a projection onto the subspace 
spanned by the {V'o-} 

/ ' Anry!t)^ '^=^^°' /^^^(r,V',i)=7^./W + ^.,(0)A,/rfv'A,/«, (D16) 
where P is a projection operatör 

VX = ug^°^a j d\ıl)aX. (D17) 



The linearized colUsion operatör of (D4) now takes the simple form 



£/(!) = {l-V) vf 



(1) 



y 



The first term represents the fact that C has a null subspace due to the moment conditions 



1 



dv 



V = 0. 



(D18) 



(D19) 



This is the usual BGK-like operatör with a single, infintely degenerate point in the spcctrum for ali functions of the 
orthogonal subspace. The second term is a projection onto a specific function in the orthogonal subspace and is the 
new efîect of the asymmetric Gaussian approximation, or the non-zero value of B,j • 
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Finally, noting that 

J dwDijDijg^''^ = ^nm''B^ (D20) 
allows the linearized operatör to be written 
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FIG. 1. Comparison of coefficients cg(q;), CGi{a) and cm(q;) with CB(a). 
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FIG. 2. Comparison of thc HCS dividcd by thc Maxwclhan for thc Gaussian model and the velocity independent collision 
frequency model with the Sonine approximation using both cg(o;) and CB{ct) for a = 0.8. 
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FIG. 3. Illustration of tlıc cross ovcr of tho HCS for the Gaussian model to the intermediate form, eq(67), by plotting 
dv* ln0(t;*) for the HCS, the asymptotic form and the Ma:x;welhan for a = 0.8. 




FIG. 4. Illustration of crossover of the HCS for the velocity independent model from the Gaussian given by Eq.(76) to an 
algebraic decay, Eq.(77) for a = 0.8. 
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FIG. 5. The exponent of algebraic decay for the Mcixwell model, Eq.(23), and velocity independent Gaussian model, Eq.(77) 
plotted as a function of a. 




FIG. 6. Plot of the fitted dimensionless viscosity rf = 77/770 for the velocity dependent and velocity independent collision 
frequency Gaussian models with the Boltzmann and Maxwell model results. 
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FIG. 7. Comparison of dimensionless thermal conductivity k* = k/ko as calculated from the velocity dependent and velocity 
independent collision frequency Gaussian models with the Boltzmann and Maxwell model results. 




FIG. 8. Comparison of n* = hu/Tko as calculated from the Gaussian models with the Boltzmann and Maxwell model 
results. 
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FIG. 9. Plot of dv* İn for the Gaussian model and velocity independent Gaussian model and the Maxwellian for a = 0.8. 
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